Late time tails in the Kerr spacetime 
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Abstract 



Outside a black hole, perturbation fields die ofi^ in time as 1/t". For spherical holes n = 2£ + 3 where 
£ is the multipole index. In the nonspherical Kerr spacetime there is no coordinate-independent 
■ meaning of "multipole," and a common sense viewpoint is to set £ to the lowest radiatiable index, 

O' although theoretical studies have led to very different claims. Numerical results, to date, have been 

controversial. Here we show that expansion for small Kerr spin parameter a leads to very definite 
fT^ ■ numerical results confirming previous theoretical analyses. 

(N 
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^ '' Perturbation fields outside a spherically symmetric black hole die ofi' in time in a way that has been well understood 
5—1 for more than 35 years This understanding is closely tied to the fact that the spherical background allows the 
I fields to be decomposed into multipoles, each of which evolves from initial data independently and can be studied 

with a relatively simple 1+1 computer code. The evolved radiation field starts with oscillations characteristic of the 
' details of the initial data, then undergoes an epoch of quasinormal ringing, and lastly falls off in time t in the form of 
, a late-time "tail" t~'^^~^, where £ is the multipole indexQ. 

QQ I For the nonspherical Kerr black hole the situation has been anything but clear since Kerr perturbations cannot be 

T-H ■ separated into independently evolving multipoles 0]. Certain symmetries do apply, however. Perturbations can be 

I separated into nonmixing azimuthal Fourier modes e*™''^ and into modes even or odd with respect to reflection through 

• the equatorial plane. This has given rise to what we shall call a "common sense" viewpoint in which approximate 
[ spherical symmetry applies to the distant radiation field, and the late-time behavior of the evolution of initial data 

• is dominated by the lowest multipole index (i.e. , the most slowly dying tail) compatible with the azimuthal and 
] equatorial symmetries of the initial data. Thus, for example, a scalar perturbation field whose initial data has m = 

^ ■ and that is symmetric with respect to the equator will, at late time, be predominately a monopole and will have the 
t'^ tail of a monopole. 

Supporting this viewpoint is the fact that, without spherical symmetry in the background, there are no preferred 
I angular coordinates. The radial r and polar 9 coordinates are mixed differently in different systems of coordinates used 
to describe the Kerr spacetime, such as Boyer-Lindquist (BL) coordinates Q or Kerr coordinates A "multipole" 
is specific to the coordinate choice, and therefore cannot determine a physical effect, like the rate of decrease of the 
field. 

Theoretical work has argued against the common sense viewpoint, and claims have appeared of numerical results 
to siipport both sides of the argument. The argument for something other than "common sense," was first given by 
Hod [a, 0, 0, 0, (see also Barack and Ori[ll|, 121). Hod considered initial data that has only a single multipole 



Yim in BL coordinates. By looking at the zero frequency limit of a Fourier transform. Hod argued that the tails of a 
massless scalar function would have the following dependence on time and on multipolarity of the initial data: 

y«m/t^^+^ i^m or £ = m+l 

j^i+m+i £ - m > 2 (even) (1) 
W mA'+"+' e-m>2 (odd). 

The common sense results seemed so compelling, that numerical work was immediately sought that would settle the 
issue, but numerical tests required rather delicate 2+1 codes. Krivan[l3| was the first to attempt this work, using a 
scalar field with an initial outgoing pulse with BL multipole indices £,m = 4,0. The Eq. ([T]) prediction for this case 
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is a monopole tail while the common sense prediction is a monopole. (This £,ni = A,0 case is the simplest 
scalar case for which there are controversial predications, and we will consider it here as the primary test case.) 
Krivan's results weakly suggested a t~^-^ law, but Krivan pointed out serious numerical problems caused by angular 



differencing. Subsequent studies of the m = 4, case, by Burko and Khanna[l4l|. and by Scheel et a /.llSll both found 



the common sense t~'^ result. But in reference lj| penetrating Kerr coordinates were used and in [15| Kerr-Schild 
coordinates were used. A constant time initial hypersurface in either of these systems is different from a hypersurface 
of constant BL time, and therefore neither of these studies represented the same problem as that to which Eq. ^ 
and Krivan's results apply. 

A new element was added to the unresolved issues when Poisson published a careful analysis [3l of tail behavior 
within linearized general relativity. Since tails appear to develop in the weakly curved background far from a central 
black hole, and to be insensitive to the details of the strongly curved inner regionpjj, it seemed reasonable that 
Poisson's analysis would be relevant to the question of tails in Kerr spacetime, at least as an approximation for tails 
that evolve from initially outoing pulses that start far from a central hole. Poisson's analysis, fundamentally different 
from Hod's, led to precisely the results in Eq. H]). 

We present here a new approach for computationally probing the late time evolution of tails in the Kerr spacetime 
and, in principle, in other nonspherical spacetimes. The advantages of this approach are: (i) it gives a clear meaning 
to "multipoles" since it uses a spherical operator for evolution; (ii) there is no angular differencing, and hence it avoids 
the errors pointed out by Krivan; (iii) the method gives convergent, clear numerical answers to the controversies of 
tails in the Kerr spacetime. The new approach expands fields and the equations that govern them in powers of the 
spin parameter a, of the Kerr metric. 

We show here the approach as applied to a scalar field ip, both for simplicity, and because previous numerical work 
has all been for a scalar field. We make a further, minor simplification by choosing the scalar field to have no 
(azimuthal) dependence. More details and the case for more general fields will be published elsewhere pj|. 

In Boyer-Lindquist coordinates, the Teukolsky equation for this case (equivalent to iJj ^ = 0) takes the explicit form 



(r^ + a^)^ 



- sin^ I 



where A has the usual meaning — 2Mr + a^, with M the mass and a = J /M the angular momentum parameter. 
Units are used in which c = G = 1. 

For any initial data, the field that evolves will depend on the spin parameter a and we expand both -0 and the 
operator of Eq. ^ in powers of a/M: 

^ = V^(o)+a/MV''^^ + (a/M)^ + L = L^") + (a/M)^ i'^) . . . (3) 

Here the ■(/'*■"'•' are functions of t, r and 9. The equations that result for the even powers of a/M are 

L(o)[V(")] = (4) 
L(o)[V'(2)] = -L(2)[V,W] (5) 
L^^i^f^)] = -LWiV/^^l-L^'^iV-^'^] (6) 

and so forth. 

The evolved fields can be solved order by order. Since the right hand sides in the above sequence are treated as 
known driving terms, only the operator L^*^^ need be inverted, but L^^'^ is just the spherically symmetric Schwarzchild 
operator, so the multipoles contained in the solutions will only be those that appear in the driving terms. 

We focus now on the primary test case € = 4, and on the question of the exponent for the late time tail. Since 
the initial data has £ — 4 there will be an evolved field that is zero order in a/M, i.e., the purely £ = A field that is 
evolved by the Schwarzschild operator i'"^. The i — A zeroth order field V'^'^^ will provide a source term on the right 
hand side of Eq. But i*^^^ (the only term in the full set of L*^"-' that is not spherically symmetric) contains a term 
-M'^sin^ed'i, and 



sin^ 9Pi{coa9) = -(10/33) P6(cos6') + (38/77) ^4(008 0) - (4/21)P2(cos 
ive a nua.drunole niece, and hence the solutio 



The driving term in Eq. ([5]) will then have a quadrupole piece, and hence i^^^K the solution of Eq. ([5]), will contain a 
quadrupole. In a similar manner, the second-order quadrupole will provide a fourth-order (in a/M) monopole source 
term in Eq. ([6|). The fourth-order field will therefore have a monopole, and hence a monopole late-time tail. 
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Paying attention only to the multipoles of interest for this late-time tail, we write the relevant multipoles of the 
relevant orders as follows: 

=r-Vf (t,r)F4(cos(?) ^^^\t,r,0) ^ f^^\t,r)P2{cos9) ij^''\t,r,6) ^ f^''\t,r)Po{cosd) (7) 
and the order-by-order 1-1-1 wave equations of interest become 

- a;./f ' + (20 + /<». = „ (8) 

- . (HM) _| If _ HM) ' . ,10, 

The numerical computations were carried out on a t, r* characteristic grid, where r* is the standard Schwarzschild 
tortoise coordinate, with At = Ar* , and with no boundary conditions. (Computations were carried out only in the 
domain of dependence of the initial spatial grid.) All three fields f^^\ ^ were evolved simultaneously. Initial 
data for jf\ at i = 0, was chosen to be a Gaussian pulse and was made (approximately) outgoing by taking the 
initial value at grid point r* to be replicated after a time step A^, at the spatial grid point r* -I- At. Initial data were 
taken to be zero for and for /q^"*, so that the final fourth-order monopole was only the result of the initial pure 
£ = 4 data. 

Typical results are shown in Fig. [TJ At earliest times these curves show the nonzero starting value of the zeroth 
order field jf'' , the field for which there is nonzero initial data. This field then "ignites" the second-order quadrupole 
which in turn ignites the fourth-order monopole. All fields go through periods of quasinormal ringing (each with a 
different frequency) around t/2M =100 and settle into the power-law tails 1/t" that show up as straight lines in the 
log-log plots. With better than 1% accuracy, we find that /^"^ has an exponent n — 11, and f^'' has n — 7, values 
that are in excellent agreement with the n = 2£ + 3 rule and that demonstrate that each tail is evolving without 
significant influence from the lower-order tail that, at earlier times, "ignited" it. 

Of particular interest, of course, is the exponent of the monopole tail, which we found to be n = 5 to high accuracy 
(better than 1%). We ran this case for several grid sizes. For Ar* less than 0.025(2M) we found some differences of as 
much as 5% in the early waveform details, but negligible difference in the computed value of n. The major source of 
error lay in determining the asymptotic value of the exponent. To do this we ran until times of about t — 4000(2Af) 
with runs that required around a day on a workstation. Exponent n was computed for segments of the late-time 
tail, and was found to be 5.0495 for the latest segment. The results showed clear evidence that the value of n was 
slowly asymptoting to a slightly smaller value. We do not consider it plausible that at sufficiently late times it would 
decrease to n = 3. 

From a numerical point of view this is truly remarkable. Since the fields at late time are evolving without source (as 
shown by the success of the 2i + 3 rule for the lower order multipoles), the monopole should have the exponent n — 3. 
Indeed, modifications in the computation do change the exponent from n = 5 to n = 3. We see this change if we put 
in nonzero initial data for /g^-* , or for /j^-* . We see this change also if we arbitrarily turn off the evolution of the f^'^'^ 
field at some intermediate time and let the /j^'' and /q^'' fields continue to evolve. The overwhelming tendency for 
the exponent to be 3 rather than 5 convinces us that there is no error we have overlooked in our program; any error 
would almost surely lead to n = 3. 

The delicacy of the the n — 5 result underscores the numerical advantages of the l-fl computations in Eqs. ([5]) 
- over 2-f 1 codes, even though the set of equations contains second time derivatives applied twice on the right 
hand source terms. The method turns out to be accurate enough that we have been able to go one step further, and 
treat the initial £,m = 6,0 case for which the common sense monopole exponent is n = 3, while Eqs. ([T]) predicts 
n = 7. We have found n — 7, again to within about 1%. 

There remains the possibility that the n = 5 result is a consequence of cancellation of contributions to the 
tail, and that these cancellations only occur to lowest order in a/M (that is, to order (a/M)'^). In a forthcoming 
paper[Tl] we shall present the results carried out to higher order. We shall also provide further details of the method 
used, and a wider set of numerical examples, including those for nonaxisymmetric initial data and for gravitational 
perturbations, and tests of Poisson's approximation [16j. 
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FIG. 1: 



Evolution of the zeroth-order 1 — 4 multipole 
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the second-order i — 2 multipole /j , and the the zeroth-order 



monopole /g for a scalar field in a Kerr spacetime. 
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